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We study the dynamics of small droplets of polydimethylsiloxane (PDMS) silicone oil on a vertical, 
perfectly-wetting, silicon wafer. Interference videomicroscopy allows us to capture the dynamics 
of these droplets. We use droplets with a volumes typically ranging from 100 to 500 nanolitres 
(viscosities from 10 to 1000 centistokes) to understand long time derivations from classical solutions. 
Past researchers used one dimensional theory to understand the typical t 1 ^ 3 scaling for the position 
of the tip of the droplet in time t. We observe this regime in experiment for intermediate times 
and discover a two-dimensional, similarity solution of the shape of the droplet. However, at long 
times our droplets start to move more slowly down the plane than the t 1//3 scaling suggests and 
we observe deviations in droplet shape from the similarity solution. We match experimental data 
with simulations to show these deviations are consistent with retarded van der Waals forcing which 
should become significant at the small heights observed. 

PACS numbers: 



I.INTRODUCTION 

The motion of a liquid drop falling due to gravity on 
a vertical solid substrate is a classic problem in fluid 
dynamics []], Q. The movement of the fluid front de- 
pends on a balance between viscous and surface tension 
forces with gravity providing the body force. In the case 
of extremely small droplets, it is expected that in the 
late stages of spreading when the drop thins appreciably, 
subtle interplay between van der Waals and surface ten- 
sion forces may become important. A detailed study of 
droplets is also technologically important in applications 
where the dynamics of droplets with solid substrates is 
used, e.x. in microfluidic devices and inkjet printing of 
organic electronic circuits |3[. 

When a uniform fluid layer (line source) flows under 
gravity, the advancing contact line develops a fingering 
instability which has been the subject of intense study 
@, H, @, 0, H, [H . In the case where surface tension can be 
neglected, the contact line shows a typical scaling rela- 
tionship with time which was analyzed by Huppert [l| for 
arbitrary angles of inclination. Since the flow of liquids is 
technologically important for various coating processes, a 
through analysis of gravity driven flow for point and line 
sources (with zero surface tension) was done by Lister [7j . 
Recently Gonzalez et al. [8] extended this study to the 
case of a finite volume fluid strips falling under the ac- 
tion of gravity. Their careful numerical and experimental 
work showed deviations from Huppert's similarity solu- 
tions due to surface tension effects. These effects lead to 
the formation of a capillary ridge in the front before the 
onset of the fingering instability. Furthermore Gonzalez 
et al. predict the wavelength of the fingering instability, 
noting it can be controlled by the surface tension. The 
shape of the profile has been investigated in detail by 
Hocking [9[, where he identifies three distinct regimes of 
the fluid sheet before the onset of the instability. 



Depending on the wetting properties of the substrate, 
the dynamics of finite volume droplets under the com- 
bined effects of surface tension and gravity can give rise 
to a number of interesting shapes 10, EH • The capillary 
number is the main parameter controlling the morpho- 
logical shape transitions and dynamics, in contrast to the 
constant flux flow of a liquid. For droplets with a finite 
contact angle, one can have corners, cusps or pearl like 
droplets depending on the capillary number. The corner 
formation in cusp-like shapes have been linked to the ex- 
istence of conical similarity solutions in the lubrication 
equation [l2j |. while recent numerical work predicts the 
existence of chaotic shedding states [Hj. 

In the present work, we concentrate on completely wet- 
ting substrate, where the contact angle goes to zero. This 
allows us to make detailed comparison between theory, 
numerical simulations and experiments, without using "a 
priori" contact angle versus disjoining pressure relation- 
ship . We capture the intermediate-time dynamics by 
a two-dimensional constant volume analysis of the lubri- 
cation equation which includes gravity and surface ten- 
sion. Only in the long time limit do the effects of the 
intermolecular forces become important. Careful numer- 
ical study of these droplets and qualitative comparison 
with experiments allows us to understand the importance 
of disjoining pressure in the dynamics [15|. We use a pre- 
cursor film in the numerical simulations [16[ and only the 
long-time dynamics is sensitive to precursor height val- 
ues. A comprehensive study of droplets under the com- 
pletely wetting scenario can reveal how modifications to 
the disjoining potential have to be made for partially wet- 
ting surfaces. It also brings out the important similarities 
and differences between a constant volume drop and a 
single finger formed by the instability of the contact line. 
Recent work on partially wetting substrates has shown 
the importance of taking these van der Waals forces into 
account for very small drops [T^ . 
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FIG. 1: Schematic of the experimental setup where rays re- 
flected from the top and bottom of the drop (inset) creates 
the interference pattern. 



In the next section, we give a detailed outline of the 
constant volume problem, along with the experimental 
setup used to investigate the problem. We will also show 
typical experimental results and how the initial circular 
shape of the drop changes from a circular plane form to 
a form with a regular scaling law and then finally to a 
nearly elliptical plane form. In Section III, we under- 
take a detailed theoretical analysis of the two dimen- 
sional lubrication approximation for drops with constant 
volume. The numerical results are compared to exper- 
imental data for position versus time, longitudinal and 
transverse shape of the drop at intermediate times. In 
Section IV, we present numerical results with van der 
Waals forcing and compare them with experimental re- 
sults when the dynamics of the front shows strong devi- 
ations from t 1 / 3 scaling law at very long times. 



II. OUTLINE OF THE CONSTANT VOLUME 
DROPLET PROBLEM 

We use Polydimethylsiloxane (PDMS) silicone oil as 
our working fluid on a pretreated silicone wafer. The 
oxidized silicon wafer is first wetted with the silicone oil 
(density 0.986 g/cc and surface tension 20.9 dyne/cm), 
then rinsed thoroughly with hexane and finally cleaned 
with a methanol wetted tissue [18]. This covers the sub- 
strate with a precursor film of PDMS which is a single 
monolayer thick and makes it a completely wetting sub- 
strate. The PDMS (United Chemical Technologies) oils 
used had viscosities ranging from 1.5 centistokes to 1000 
centistokes. Since the timescale for the various effects are 
controlled mainly by the viscosity, changing the viscosity 
allows one to follow the change of shape from a circular 
to an elliptic cross section. The volume of the drop was 
varied from around 100 nanolitres to a few microlitres by 
using a micropippette (Biohit Proline Micropippetors). 

The droplet was placed on the surface of a silicon wafer 



FIG. 2: Close up of Interference fringes from PDMS oil 
droplet of viscosity 100 centistokes and volume 0.85 mi- 
crolitres. Note the Capillary ridge at the leading edge. 



(which was mounted on a vertical plane) and was illu- 
minated by a collimated He-Ne laser beam (wavelength 
632.8 nm) as shown in Figure [TJ The rays reflected from 
the top surface of the oil drop and the surface of the wafer 
gave a typical interference pattern as shown in Figure [2j 
The position and shape of the front as well as the fringes 
can be followed as a function of time. It should be noted 
that most of the analysis will be concerned with the main 
body of the drop. Surface tension effects create a capil- 
lary ridge at the tip of the drop. The back of the droplet 
moves extremely slowly on timescales compared to the 
front of the drop. 

From the images (Figure [2|) three main sections of the 
drop can be distinguished, a dark front which comprises 
the capillary ridge, the main body of the drop where most 
of the interference fringes can be seen and the rear section 
of the drop which moves extremely slowly on the interme- 
diate and long time scales. In addition to following the 
front of the drop, the interference fringes (Figure [2]) will 
allow us to draw some qualitative conclusions regarding 
the surface profile, which can be compared with numeri- 
cal simulations and theory. 

We note that the height of the fluid is much smaller 
than the extent of fluid in the plane, and the slope of 
the free surface is small except near the capillary tip. 
Hence we can use the lubrication approximation for an- 
alyzing the motion of the bulk of the fluid. It is also well 
known that if the Bond number, B = pgW 2 /a, where 
W is a typical transverse length scale, is much greater 
than 1, then the effects of surface tension and contact 
line motion can be neglected. Depending on the the vis- 
cosity and surface tension our typical Bond numbers are 
around B ~ 0.2 — 0.8. Hence we are in the regime which 
are in the opposite regime of those analyzed by Hup- 
pert and Lister [I], 0]. In the case of a droplet on a 
partially wetting substrate a function of the Bond num- 
ber also appears as a control parameter which determines 
the transition between shedding and nonshedding states 
[l3( | . The capillary lengths for our PDMS oils are around 
\J a/ pg ~ 0.1 — 0.2mm where a is the surface tension and 
p is the density of the drop. Assuming a hemispherical 
drop of volume around 0.2 — 0.5 microlitres, the initial 
radius about the same scale as the capillary length. We 
choose the free scaling parameter, the height scale, to be 
this value. 
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When the drop is initially put on the substrate, high 
curvature at the contact line pushes the fluid radially 
outwards. Depending on the viscosity of the drop the 
circular plan form maintains itself until the gravitational 
forces overcomes inert ial and capillary effects. If the ini- 
tial radius of the drop is R and the typical length scale 
in the downslope direction is given by L, then the sur- 
face tension driven spreading, which is dependent on the 
curvature k, is given by k ~ R/L 2 . As was shown in 
recent experiments of Biance et al. [HI, m the inertial- 
capillary regime the length should scale as ~ VDt, where 
the dynamic diffusion coefficient is given by D = y/crR/p. 
We do not observe any y/t kind of dependence in the 
timescales when we start taking the first measurements 
(1 — 2 seconds). This is borne out by the fact that the 
dimensionless Ohnesorge number, Oh = rj/^apR, which 
distinguishes between the inertial and the viscous regime 
[2Q( | . is around 1 for the viscosities rj in our experiments. 

In Figure [3j we show typical examples from the three 
main regimes we see in the experiments. Initially, the sur- 
face tension forces dominate and we have an nearly cir- 
cular drop. In the case of a purely gravitationally driven 
flow, neglecting the effect of surface tension the initial 
axisymmetric spreading was found to be proportional to 
t 1 ! 2 by Lister [7[ and gave away to predominantly downs- 
lope flow when t ~ 0(1). 

In the second regime we see a very regular motion 
where both gravity and surface tension have to be taken 
into account. The transients typically die out in the time 
scale « y/rj. Looking at the position of the tip as it varies 
in time, we see t 1 ^ 3 scaling law relationship (Fig. [6]) for 
the intermediate times. Huppert [1] first explained this 
scaling law using one-dimensional theory without surface 
tension. We extend this theory to two dimensions and 
use asymptotics to allow a detailed understanding of the 
transverse features of the droplet. 

At very long times, we start to see deviations from the 
intermediate-time scaling laws. The motion of the tip 
slows down significantly, and the droplet starts to spread 
laterally at an increasing rate. We note through other 
works |Tl|, E3 that at very small height ranges near the 
end of our runs, van der Waals forces should start to be- 
come important. Looking at the effects of these forces 
in our numerics we can see the same shape changes as 
we see in the experiments. In particular, we can match 
the deviation in tip position from the t 1 ! 3 scaling law in 
experiment and simulations. From this evidence we con- 
clude that van der Waals forces are the cause of long-time 
deviations from the classical solutions seen at intermedi- 
ate times. 

Ill FORMULATION 

We study the gravity-driven lubrication flow of a small 
droplet of Newtonian fluid on a vertical solid flat plate. 




FIG. 3: The three regimes of droplet motion shown in ex- 
perimental footprints and corresponding numerical profiles. 
The first pictures are from the initial transients. The second 
set corresponds to the important asymptotic regime where 
the gravitational effects are balanced by inertial and surface 
tension effects. Finally, as the height scale of the drop be- 
comes much smaller we enter a third regime where Van der 
Waals forces become import ant. See section (III) for detailed 
explanation of numerics. All the numerical heights are in di- 
mensionless units. 



The finite-mass droplet has a constant density p, viscos- 
ity p and surface tension a, and the flow is driven by 
gravitational stresses and interacts with the fully wet- 
ting surface. We define a coordinate system with the in- 
coordinate in the direction of gravity and the z-direction 
normal to the plate, as shown in Figure [H The y- 
direction is transverse to the direction of flow and com- 
pletes the right-handed coordinate system. We denote 
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FIG. 4: Geometry of the problem 

the height of the free-surface of the flow as z = h(x, y, t). 

We shall now briefly review the derivation of the lubri- 
cation model including long range molecular forces [211 ]. 
We rescale the dimensional (hatted) lengths, time, ve- 
locity, and pressure to yield corresponding dimensionless 
quantities, 



x = L x. y = hy 1 z = Hz, t=—t, 



u = XJu, v = XJv, 
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and use these rescalings on the the Navier-Stokes equa- 
tions for the viscous, incompressible flow where the modi- 
fied pressure P directly incorporates the hydrostatic pres- 
sure and the van der Waals molecular interactions, 



P=p + xSt~ 1 +Han(ft). 



(1) 



The dimensionless parameters present are the aspect ra- 
tio, the Stokes and Hamaker constants as well as the 
reduced capillary number, 
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where g is the acceleration due to gravity and A! is the 
dimensional Hamaker constant. 

At z = 0, the no-slip boundary conditions, u = w = 0, 
are applied. At the free-surface, z = h(x,t), the bound- 
ary conditions are the kinematic stress balance condi- 
tion, the normal stress balance and the tangential stress 
balance. We follow the derivation in Oron, Davis and 
Bankoff [2l| integrating the Navier-Stokes equations to 
find the depth averaged velocity and then using the ki- 
netic boundary condition to produce a partial differential 
equation for the height, 

= HaV • (h 3 \/U(h)) - • (h 3 VV 2 h). (2) 
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where V = (J^, J^) is the gradient on the plane. 



Despite the relatively small constant in front of the van 
der Waals term in (|5J), we must keep this term to under- 
stand the longtime dynamics, as we will show in section 
IV. We will use l/h 4 form for the van der Waals force 
n(/i), which comes from the dielectric properties present 
in the liquid-solid layered system when the liquid has 
no ionic species [22|. We choose to incorporate the as- 
pect ratio into the Hamaker constant because the van 
der Waals interactions are important near the contact 
line at the edge of the droplets. This term then must be 
included to correctly model the physics despite the small 
value of the Hamaker constant in the lubrication limit 
(e — > 0). As the film thins the 'retarded' or long-range 
part of the van der Waals interaction becomes impor- 
tant. Casimir and Polder first showed that in the 
limit of long distances the interaction energy for the van 
der Waals interaction changes from 1/r 6 to 1/ r 7 .In order 
to calculate the disjoining pressure from this functional 
form one has to integrate the two parallel domains using 
the full Lifshitz theory, which has been done by Derjauin 
et al. [22[. This 'functional' form of the van der Waals 
interaction comes into effect when the film thickness is 
above 10 nanometres. The retarded van der Waals inter- 
action has been used in modelling the dynamics of thin 
tear films in the eye amongst others [24|, [25(. For the 
disjoining pressure of the form A' /h 4 , A' has dimensions 
of erg-cm, whereas for the non-retarded form, A' /h 3 , A 
has dimensions of ergs. 

We now choose our scaling similar to [26], 



99 ) 



1/3 



3/iL 



(3) 



which sets all the coefficients in front of the terms in ([2]) 
to unity except for the Van der Waals interaction term, 



dh dh 3 A ^ /o 1 



V-(/i 3 VV 2 /i). (4) 



We keep the ability to adjust or change the strength 
of the van der Waals interactions though the non- 
dimensional parameter, 
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(5) 



so we can understand the effect of the intermolecular 
forces on the fluid. Our other major parameter in the 
problem will be the volume of the droplet V which when 
scaled becomes V = 7777- However, for our comparisons 
between numerics and experiments we can avoid deal- 
ing with this parameter by setting our free height scale 

T 2 

H = ^7, and therefore obtaining a scaled volume of ~ 1. 

One-dimensional numerical calculations were pre- 
formed using standard Crank-Nicolson framework [271 ]. 
All two-dimensional numerical solutions for the govern- 
ing equation (|4|) and the scaled version of this equation 
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we discuss later ([TT|) u se a midpoint, alternating direction 
scheme as seen in [28]. These methods along with New- 
ton's iterations for the non-linearity, allow for second- 
order accuracy in time and space while keeping reason- 
able execution time. 

In all cases, we must estimate the pre- wetting height of 
film used in the experiments as this value is not directly 
measurable. This was chosen to be two to three orders of 
magnitude less then the smallest heights achieved by the 
droplet. Throughout this range we find the same quali- 
tative results seen in the next sections for scaling powers 
and droplet shapes before the van der Waals forces be- 
come important. As the film height becomes extremely 
small, however, the strength of the van der Waals term is 
affected by the wetting height. This can be understood 
by looking at the terms in the governing equation (j4]), 
and noticing that most of the terms depend on positive 
powers of h while the van der Waals has a negative power 
of h and therefore is more sensitive to small changes in 
h. Pre- wetting dependence for the quantitative motion 
of fluid has been hypothesized by many authors 0, @, [2l[ 
as explanations for large deviation in proportionality con- 
stants while still having very good fits for scaling powers. 
We clarify the cause of these deviations in the Appendix. 



valid from an arbitrarily chosen x = to the tip at x t - ip = 
(^Y^-Va) 2 / 3 ^/ 3 determined by conserving the profile area 
Va- The solution is zero for x outside this region. The 
appearance of a "shock" or "jump" is not too surprising as 
we have removed the fourth-order smoothing. However, 
the main effect of this smoothing is to produce a capillary 
ridge, which appears due to the finite surface tension of 
the fluid, that smooths out the shock at xu p . 

The numerics show the predicted square root scaling 
(|8j) of the profile height (Fig. [5]), which can be compared 
with the experimental data from analyzing the spacings 
in interference fringes. This shape of the the profile is 
markedly different from the profile in the case of gravity 
driven fingers from a line source [26[, or in the case of 
Marangoni driven fingering [18[- For Marangoni driven 
fingers or for gravitational fingering from a line source, 
one has a very flat intermediate section followed by a 
sharp tip at the end. 
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A One-Dimensional Solutions without van der Waals 
Forcing 



We start the analysis of one dimensional solutions 
without van der Waals forcing by taking the governing 
equation (|4]) in one-dimension, for solutions that are uni- 
form in the ^/-direction and with no wetting i.e. A = 
we obtain the governing equation 



dh dh 3 _ 3 

dt + dx ~ [ xxx)x ' 



(6) 



This equation has been studied extensively for both con- 
stant volume [l|, [|[ and constant flux [7] cases and we 
will review the main results here. We verify our results 
in intermediate times before van der Waals forces be- 
come important using theoretical results from this equa- 
tion. Also, we use the scaling laws presented to under- 
stand the novel, two-dimensional asymptotic solution in 
the next section. 

Ignoring the fourth order smoothing for the moment, 
we obtain the same equation that Huppert first used [1] 
to describe the motion of a droplet, 



0. 



with the solution 



dh dh 3 
dt dx 



h A =o(x,t) = x j — 



(7) 



(8) 
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FIG. 5: A numerical profile h(x,y = 0) fit with Cx 1 ^ 2 and 
the analysis of the interference fringes from a 1000 centistokes 
run. Fringe profile is fit to a quadratic function in the inter- 
mediate section of the drop away from the tip (shown for the 
first 25 fringes). 



The scaling law for the tip position x t [ p oc t 1 ! 3 has been 
seen throughout the literature for gravitationally driven 
fluids [3, Q. In our experiments and numerics in one 
and two dimensions, we see this scaling for intermediate 
times. In Figure El we show typical experimental results 
for 50 and 100 centistokes of assorted volumes along with 
a numerical run with A = 0,to show the regime of t 1 ! 3 
behaviour . The numerical run has, as expected, a nearly 
perfect t 1 / 3 scaling, however, we see deviations from this 



6 



scaling in the experiments at long times which will be 
attributed to van der Waals forces in later sections. 
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FIG. 6: Typical data sets for viscosities of 100 or 50 centis- 
tokes and volume ranging from 100 nanolitres to 1.0 micro- 
liters graphed with a reference numerical solution (solid) for 
A — showing the t 1//3 scaling. Squares and diamonds corre- 
spond to 50 centistokes, others are 100 centistokes. Note the 
t 1 / 3 scaling in the data at early times and long time deviations 
from this scaling in all cases. 



B Two-Dimensional Solutions without van der 
Waals Forces 

In the Appendix, we derive a two-dimensional, long- 
time, asymptotic (near the tip), series solution to our 
governing equation (j4]) 

«*.».«> -^(l-SjSjO+tV) <» 

where w(x, t) = W(l - z£" 1/3 /£ti P ) 1/4 is the width of the 
droplet. Note, our solution still has two constants that 
need to be fit W and £tip, due to the lack of a matching 
asymptotic solution for the back of the drop (near x = 0). 
The first constant W controls how quickly the footprint 
width scales with the distance from the tip. The sec- 
ond constant £ t i p , in the physical variables, becomes the 
scaling constant for the tip position x = £tip^ 8 ^ 15 ^ 3 - 
Numerically, we have noticed that these constants are 
quite independent of initial droplet shape as long as the 
initial volume is contained within an area that is about 
half of the eventual width W. However, as we look at 
the experimental values for £ t i P from Figure [6] and W 
from footprint data in Figure We find quite a bit of 
spread in these constants despite the very close agree- 
ment for the scaling. This spread in the values has 
been seen in other droplet literature [H, Q and is dis- 
cussed in the Appendix. The asymptotic solution is not 
valid near the tip of the droplet. The center line height 
(y = 0) matches the square root scaling discussed earlier 



in the one-dimensional theory ([8]). Also, if we look for 
the point where the width of the drop w(x,t) becomes 
zero (the tip of the droplet), we see clearly the tip scaling 
x t ip oc t 1//3 , as discussed in the previous section. In ad- 
dition to matching the expected one-dimensional theory, 
we have a novel explanation for a number of the trans- 
verse features we see in experiment. 

We clearly see the inverse parabolic cross-sections pre- 
dicted by © (for constant x and t) in the numerics as 
well as the experiments (Fig. [7]), when data is taken away 
from the tips of the droplet. In the Appendix, we show 
that the 0{y A ) term in the series solution must be zero 
therefore our parabolic fits are nearly perfect. Experi- 
mentally, these parabolic cross-sections are seen for in- 
termediate times before the van der Waals term becomes 
important. The final important feature we can predict 




FIG. 7: A typical numerical cross-section (taken at an x value 
of about 80% of the total drop length) for a solution to j4j) 
with A — (top) and a typical experimental cross-section 
(taken at a postion about 80% of the total drop length) at 
intermediate times fit to a parabola 



using the asymptotic solution (O, is the shape of the 
footprint of the droplet (where the height becomes zero). 
The theory predicts a l/4th scaling law for the width w of 
the droplet as a function of position along the droplet x. 
Our numerical solutions for A = show the correct 1 / 4 th 
scaling in the width of the droplet in the footprint, esti- 
mated here by a contour of height h = 0.05 (Fig. [8]). Note 
that we see some deviation at the tips of the droplets. 
The deviation for x near zero is not surprising as we 
derived our solution asymptotically near the tip of the 
drop. However, we see that the solution is surprisingly 
valid through most of the droplet length which again can 
probably be attributed to the zero 0(y A ) term. Also, in 
the appendix we have scaled out the x derivatives near 
the tip which should become relatively small over time. 
However, as our numerics are not steady st ate solutions, 
but merely long time solutions, and as the x-derivatives 
remain strong right near the edge of the drop, we see the 
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l/4th power law cuts off just before the tip. 
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FIG. 8: The top half (77 > 0) of a numerical footprint from a 
solution of J4} under scaling ([TO]) fit with a scaling constant 
of l/4th (top) and an representative experimental footprint 
(volume 0.1/i/, viscosity 100 cSt) fit to a scaling function. As 
with all of the experimental runs for an intermediate times 
where van der Waals forces are not yet important. The scaling 
constant (in this case b — 0.23) was close to the predicted 
value (1/4). 



IV THE EFFECTS OF VAN DER WAALS 
FORCING 

In this section we will analyze the effects of the van 
der Waals term on the shape of the droplet. The relative 
strengths and applications of the 'retarded' and 'non- 
retarded' Hamakar constant has been discussed in detail 
by French (29[. When the completely wetting droplet 
thins, the tip can typically reach a height of less than 
100 nm, which is typically in the range where retarded 
Hamaker constant will be effective. In fact, for horizontal 
droplets ellipsometric measurements during the late stage 
of spreading show an advancing 'precursor' layer which 
is about 30 nm in thickness [16[. In recent experiments 
on spreading of PDMS on partially wetting substrates, 
for drop volumes of about 100 nanolitres there was a 
crossover from surface tension to van der Waals spreading 
[T3|. Extensive work on crossover from 'non-retarded' to 
'retarded' forces for liquids on quartz systems have been 
analyzed by Teletzke et al. [lo[ . 

To understand the effects of the van der Waals forc- 
ing we turn to a change of variables suggested by the 
scaling in the tip position x ~ t 1 / 3 . Since our motion 
is mostly one-dimensional at long times, we choose to 



scale our spatial variable £ to keep the tip position con- 
stant. This change of variables is a modified version of 
a scaling proposed by Hocking [3l|, from which we can 
numerically find a similarity solution to our problem (see 
Appendix) without the van der Waals forcing. By mov- 
ing into the scaling variables from the similarity solution 
we can easily observe the changes caused by the van der 
Waals term We then pick the rest of our variables, rj the 
transverse scale, s as our "long" time variable and <fi the 
scaled height, to asymptotically balance the governing 
equation in t. 

i = t-^x 



(/>(£, 77, s)=t^ 3 h(x,y,t) 

We apply this change of variables to the full governing 
equation with wetting < i < 1 and obtain the fol- 
lowing partial differential equation for the scaled height 

- (4> 3 <Pvvv)v - 0(e~ 2s / 3 ). (11) 

For very small A our solution will evolve into to 
the two-dimensional solution © for intermediate times 
(Aexp(4s/3) = At 4/S <C 1), but it is clear from the 
scaling that the van der Waals forces will eventually be- 
come dominant. Importantly, the van der Waals forcing 
in the transverse direction becomes important much ear- 
lier t ~ A -3 / 4 while the x-direction term kicks in at 
t « A -3 / 2 . From this scaling difference we expect to ob- 
serve transverse spreading in the experiments. Note the 
numerical footprints in Figure [9] spread to approximately 
the footprint of solution with A = but then continue 
to spread in the transverse direction over longer times. 
We see the same qualitative changes in the experimental 
footprint (Fig. [T0|) as well as a loss of the l/4th scaling 
of the width predicted in the Appendix. 

The lateral spreading causes deviations in the shape 
of the drop from the similarity solution without van der 
Waals forces ((9]). Through numerical solutions we ob- 
serve the transverse spreading slightly reduces the height 
of the profile height from the values predicted by the 
similarity solution. We also see the cross-sections are 
no longer parabolic to \0(y 6 )\. These deviations are 
too small to observe in the experiments. However, we 
can easily observe a reduction in the speed at which 
the tip of the drop moves down the plane in the un- 
sealed variables as shown in (Fig. [TTj). These deviations 
from the similarity solution start at the correct time scale 
Aexp(4s/3) = At 4 / 3 « 1. We can fit various data sets 
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FIG. 9: Scaled numerical footprints of the droplet at various 
long times for A = 4 x 10 -12 . Note the solution starts near 
the A — similarity solution (solid) at intermediate times 
and we can see significant shape change at increasingly long 
times (in the direction of the arrows) due to van der Waals 
forces. 




10 




10000 100000 





















10000 




100000 



FIG. 11: A data set with viscosity 50 centistokes 
(above) graphed with the closest fit numerical solution A — 
2 x 10~ 12 ± 1 X 10~ 12 and wetting height 6 x 10" 5 .A similar 
data set with viscosity 100 centistokes (below). 
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FIG. 10: Scaled experimental footprints (100 centistokes) at 
an intermediate time (1550 seconds, diamonds) and a rela- 
tively late time (14500 seconds, crosses) showing qualitatively 
the same transverse spread seen in the above numerical runs. 



and compared their deviation in the tip position from the 
scaling law x t ~ t 1 ^ 3 . Using numerics with various val- 
ues of A we can find an approximate strength of the van 
der Waals forcing for that experimental run as shown in 
Figure ITTi 

However, in obtaining this value for the dimensionless 
Hamaker constant A we had to make many assumptions 
that could significantly effect the final value. Some com- 
plications arise in fitting stem from the undetermined 
constants in the solution (j9j). As initial lateral spreading 
and mass distribution can vary within the experiment, 
we choose to fit the data using an effective volume in- 
stead of the measured volume. The effects of the wet- 
ting height (which is unmeasurable), are even more sig- 
nificant as small changes in the wetting height can have 
fairly large effects on the deviation seen and therefore the 
matching value for A. Also, the deviation is similar for 
different functional forms for the van der Waals pressure 
n(/i). For the forms of the pressure tested in the numer- 
ics (11(h) = 1/h 3 , U(h) = 1/h 4 , 11(h) = d/h 3 + C 2 /h A ) 
we see very similar qualitative results for the deviation 
from the l/3rd scaling law. In particular, we find that 



the van der Waals forces are the contributing factor for 
the long time deviations from the expected t 1 ^ 3 scaling 
law for the tip position. 



V SUMMARY AND CONCLUSION 

As noted earlier, the scaling exponents fit very well 
for intermediate times, however, in our experiments as 
well as the experiments of others the constants associ- 
ated with these scaling laws vary rather widely. Other 
authors 0, H, HH have hypothesized that wetting effects 
contribute to this spread in values. From (fT9j) in the Ap- 
pendix, we can see a possible method for wetting effects 
to cause this spread. A matching asymptotic solution 
near the back of the drop (£ = 0) would almost cer- 
tainly depend on the form and strength of the wetting, 
as the height is very small even at intermediate times. 
This matched solution would determine the value of W 
and therefore (from conservation of mass), the £ t ip- So, 
the van der Waals form could well affect the scaling con- 
stants strongly even though they don't explicitly affect 
the droplet shape for intermediate times. 

The experimental and numerical evidence clearly 
points to van der Waals forces as the cause of long time 
deviations from classical solutions for the motion and 
shape of microliter-sized, liquid droplets. While we derive 
a asymptotic, series solution for the shape of a droplet 
without van der Waals forces ([9]) , the solution could use a 
matching asymptotic solution near the back of the drop. 
This would likely give insight into into the spread in data 
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seen in this paper and others on droplet motion. Also, 
asymptotics for small amounts of van der Waals forc- 
ing on this complete solution could complete the under- 
standing of the long-time droplet shape changes. Mea- 
surements of the precursor height, while not possible with 
current apparatus, could well become possible in the near 
future. This would perhaps allow close fitting of the di- 
mensionless Hamaker constant and test the correct form 
for the van der Waals pressure term. 



APPENDIX: APPROXIMATE SOLUTION FOR 
THE DROPLET WITHOUT VAN DER WAALS 
FORCES 

To understand the shape of the droplet we work in two 
dimensions, we turn to our change of variables ([T0|) . Our 
governing equation (|4]) under this rescaling becomes 

^ s + (<P 3 - £0/3) € = -(^^ vvv ) v + 0(e" 2s / 3 ). (12) 

Long-time numerical simulations to this equation, show a 
single long time steady solution with gravitational forces 
(j) 3 balancing the term produced by the shrinking of £ 
scale in time (Fig. |4] moved back to standard variables). 
This solution is independent of initial condition shape 
as long as the initial droplet is reasonably concentrated 
(the initial width is less then half of the final transverse 
width). This condition was not a concern as the initial 
droplet sizes in the experiments were significantly smaller 
than the final droplet in width. 

We start looking for an analytic, long-time form by re- 
moving the time dependence <p s = 0. For our analysis, 
we call £ = £tip to be the steady-state, tip position which 
will end up being close but not quite the actual tip posi- 
tion as we will ignore the capillary ridge seen in Figure O 
We will define a new similarity position variable £ as the 
scaled distance up the drop from this point £ = £ t ip- 

5tip 

In a study of fingering under gravity and surface ten- 
sion gradients, Witelski and Jayaraman (private com- 
munication) introduce a mixed-spatial similarity variable 
they use to understand a balanced fourth-order smooth- 
ing perpendicular to their driving forces. Since we also 
have a driving force perpendicular to our smoothing we 
use the same similarity variable fj to modify our govern- 
ing equation 



T] 



(i-£/6ip) a e*" 



(13) 



The similarity constant a takes a particular importance 
as the factor determining how the width of the droplet 
scales as a function of the distance from the tip (Fig. EJ). 



These new variables when put into the steady version 
of ([T2]) with A = produce the following equation for 



(i - l)h 



<P (14) 



which appears quite daunting until we realize that we 
are interested in small ( C 1 near the tip of the drop. In 
which case, we can set our scaling constant a = 1/4 and 
look at the order 0(l/£) equation 



^777777)77 H~ 



1 



4£ti P 



12 ' 







(15) 



To find an approximate solution to this equation, we 
note that the problem is symmetric about the line 77 = 
and look for a series solution to the problem 



<K£, fj) = + W 2 f + W A f) A + ...) 



(16) 



from which we find that after collecting terms of order 
ff that W4 must be 0. The full series solution for this 
equation does not seem to have an closed-form explicit 
solution. However, since the equation skips three orders 
between the non-zero W2 and wq, we can hope to get 
a very reasonably accurate solution using only the first 
two terms of the expansion (W4 = wq = wg = . . . = 0). 
Plugging in this shortened expansion into the full partial 
differential equation ([T4|) and collecting terms of order 
ff we obtain just the terms from right hand side of that 
equation 



f(0 



(17) 



This equation can be integrated once in £, but we must 
set the integration constant to zero to keep the solution 
from becoming infinite as £ — » 0. At this point we can 
just take the real solution of the cubic 



m = 



&p(i - 



(18) 



and find the form for the solution along the center line 
77 = 0. 

To get a value for W2 we would have to match this 
solution to an asymptotic solution near the base of the 
droplet (where £ <C 1) or perhaps multiple asymptotic 
solutions. As it stands, efforts to produce such asymp- 
totic solutions for the remainder of the droplet have not 
been successful. However, we can leave W2 as a fitting 
parameter like £tip an d fit solutions to the numerics and 
the experiments. To make these comparisons, we move 
back into the scaled variables 



<t>(H,V,A = 0) 



ur ' 



(19) 
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where w = W(l — ^/^tip) 1 ^ 4 and we have assumed from 
the numerics W2 is negative and for convenience changed 
our fitting constant to W = y^i^l- This series solution 
of [T2l (to 0(y 2 )) is defined for -w(£) < rj < w(£) and 
< £ < ^tip and the complete solution is defined to be 
zero elsewhere. When this similarity solution is scaled 
back to the h(x,y,t) variables using (jTUJ) we obtain the 
two-dimensional solution ([9]) and its associated equation 
for the width of the droplet. 
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